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We discuss the renormalization of the initial value problem in Nonequilibrium Quantum Field 
Theory within a simple, yet instructive, example and show how to obtain a renormalized time 
evolution for the two-point functions of a scalar field and its conjugate momentum at all times. The 
scheme we propose is applicable to systems that are initially far from equilibrium and compatible 
with non-secular approximation schemes which capture thermalization. It is based on Kadanoff- 
Baym equations for non-Gaussian initial states, complemented by usual vacuum counterterms. We 
explicitly demonstrate how various cutoff-dependent effects peculiar to nonequilibrium systems, 
including time-dependent divergences or initial-time singularities, are avoided by taking an initial 
non-Gaussian three-point vacuum correlation into account. 


INTRODUCTION 


Quantum Field Theory out of equilibrium has received 
a lot of attention in recent years, especially within the 
framework of the Kadanoff-Baym equations [1] . Coupled 
to a resummation such as the one provided by the two- 
particle irreducible effective action [5], these equations 
elude the secularity problem and yield a controlled time 
evolution even at late times [am]. Despite this major 
progress and the related growing number of applications 
including condensed matter, early universe cosmology 
or heavy-ion collisions, a consistent formulation of the 
initial value problem in the case of (super)renormalizable 
theories, that deals with the elimination of ultraviolet di¬ 
vergences at all times, is still to be constructed. Many 
interesting approaches have been devoted to understand¬ 
ing and tackling the problem, from studies in the Hartree 
approximation or in perturbation theory [SHini which 
however do not capture thermalization or are not free 
of secular terms, to approaches based on appropriately 
chosen external sources mm or on the use of informa¬ 
tion about the time evolution prior to the initialization 
time |13] which depart however in spirit from the strict 
initial value problem, and restrict the control over the 
initial state. 

In this letter, we present a consistent formulation deal¬ 
ing with both secular terms and UV divergences, in which 
the only ingredient is a proper description of the initial 
state. We focus on a simple setup that allows us to ex¬ 
hibit the features which render renormalization out of 
equilibrium difficult, while admitting an analytical treat¬ 
ment of divergences. Our main result is that nonequi¬ 
librium initial states encompassing (a particular subset 
of) non-Gaussian vacuum correlations, together with the 
usual vacuum counterterms, ensure a manifestly finite 
evolution. Furthermore, we find that initial correlations 
play a role in the elimination of divergences across all 
time-scales. 


NONEQUILIBRIUM EVOLUTION EQUATIONS 

We consider a theory involving two real scalar fields (f 
and X in d -b 1 dimensions with trilinear coupling, 

^ - ^rnl(p^ -b - \rnlx^ - ^Lpx^ . (1) 

We are interested in obtaining a properly renormalized 
non-equilibrium evolution for the correlators 

F{x,y) = ]^{{T{x),ip{y)}), p{x,y) = i{[(p{x), ip{y)]). 

( 2 ) 

While our findings can be generalized, we assume for sim¬ 
plicity that X is kept close to equilibrium by some fur¬ 
ther (unspecified) interactions, and analyze the super¬ 
renormalizable case d = 4, which admits a non-trivial 
continuum limit. 

Starting from an initial state at time t = 0 described 
by a density matrix p, the time-evolution of any ob¬ 
servable O can be obtained from the closed-time path 
representation of (O) = Tr(pO). A general density 
matrix can be parameterized by its matrix elements in 
the basis of eigenstates of the field operators at t = 0, 
(t)°-{{),x)\(t)) = 4>^{x)\(j)), where ())“ = p,x, 

{(j>+\p\(j)_) = exp (j.'^ ^ f d‘^Xi---d'^Xn 

n>0 ai,ei—±: 

• ■ ■ , X„){xi) ■ ■ ■ (l)'^l^{Xn)y (3) 

Terms with n < 2 encode the initial conditions for the 
one- and two-point functions. Non-Gaussian initial cor¬ 
relations (n > 3) enter the diagrammatic expansion of 
the generating functional Z = J 'D(j){(j)^\p\4>_)e'^ 
and can be viewed as effective n-point vertices which are 
supported at t = 0±, the upper and lower boundary of 
the closed time path C, respectively [T4|. The full self¬ 
energy for p can be decomposed as [14] 

n{x,y) = UF{x,y)-yignc{x°-y°)Up{x,y) 
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FIG. 1: One-loop contributions to the ifi self-energy Yl{x,y). 
The box represents the initial three-point correlation iaz. The 
first and second diagrams contribute to ^f/p and 11;^“^ re¬ 
spectively, while the other two are oc 5s/a(xo) and not needed 
in what follows. 


where sgn(j(x*^ ~y^) is the signum function on the closed¬ 
time path, Ss/ait) = iS{t — 0^) ± d{t — 0_))/2 and we 
omitted terms oc Sg/ai^^) which will not be needed. The 
self-energies vanish for Gaussian initial states, and 
contain an a-vertex attached to the right leg in a dia¬ 
grammatic expansion (c/. Fig.[^. For a spatially ho¬ 
mogeneous state, it is convenient to use the momentum 
representation = f F{x,y). 

The time-evolution starting from a general initial state 
at t = 0 can be described by the Kadanoff-Baym equa¬ 
tions umi] 

[d'i + vl+ 5ml)Fp{t, t') = f dt"IlF,pi^, t")pp(t'', t') 

Jo 

+n^Fl{t)Fp{o,t') 

+in^“(t)Pp(o,t'), (5) 

(d? + ^p+ Sml)pp{t, t') = dt"Ilp^p{t, t")pp(t'\ t '), 

with Vp = + TO^. The initial state enters via the 

initial conditions Fp(0,0), dtFp{0,0) = dt'Fp{0,0) and 
dtdt'Fp{0,0) as well as the non-Gaussian initial corre¬ 
lations a. The corresponding initial conditions for pp 
are fixed by the Equal Time Commutation Relations 
(ETCR) to be Pp|t=t' = dtdt'Pp\t=t' = 0 and dtPp\t=t' = 
—dt’Pp\t=t' = 1, and we assume vanishing field expecta¬ 
tion values. 

We take an initial three-point correlation into account, 
and approximate the self-energies at one-loop with cutoff 
A (c/. Fig.[l| ). The one-loop expressions for fl^^p are well 
known |15j We have in addition 

l^^{DF,,imDF,p-,it,0)zar 

-\Dp,,{t,0)Dp,p_,{t,0)zaD , ( 6 ) 



with = 1, = ±1 transformed into spatial Fourier 

modes. 

As mentioned above, we assume that the field x is 
kept close to equilibrium at all times, with correlators 
givenhy Dp, = {1 + 2nq) cos{ojq{t — t')) / {2ojq) and 
Dp,q{t, t') = sva.{ujqit — t'))/ujq, where Uq = — 1). 

This corresponds to neglecting the backreaction of the 
thermal bath to which the field p is coupled. In this ap¬ 
proximation, the spectral function is given by the equi¬ 
librium solution pp{t,t') = pp{t—t'), while the statistical 
propagator Fp{t,t') approaches the thermal solution for 
late times t, F —)■ oo m- 


EQUILIBRIUM PROPERTIES 


Before discussing the renormalization of Fp{t,t'), we 
briefly discuss the renormalization of the spectral func¬ 
tion by usual vacuum counterterms. We focus for simplic¬ 
ity on Pp'pQ for = 0, but the relevant UV properties 
extend to non-zero p,T or m^. The Fourier transform 
/Op=o(i^) = / is obtained from the vacuum 

Euclidean propagator GyadiPb) = 1/(^5 + -b Snip -b 
^vaciiPb)) as Pp‘Fq{uj) = -2ilmGyacdP5 ^ ui + ie). At 
one-loop, we find 


^vaci^'^PO 


647r2 


A - ^Arctan f) 
2 \P5) 


(7) 


There is a linear divergence which is absorbed by a mass 
counterterm = (A^/647r^)A. Then the Euclidean 

propagator admits the continuum limit G'^dvd — 
l/(T5+”^^+7b5|)) where we introduced 7 = A^/(2567r). 
The spectral function in the continuum limit reads 




—2ixu} 

(o;^ — -b 7^0;^ ’ 


( 8 ) 


and obeys J ^wPp“q(w) = —i in agreement with the 
ETCR. We note also that the spectral function behaves 
like l/w^ at large |a;|, a property that we shall use in the 
next section. This property extends to T > 0 since the 
thermal contribution to the one-loop self-energy behaves 
like }?T^jud at large external Matsubara frequency a;„. 

For the discussion below, it is finally important to real¬ 
ize that even though Gyadwd) admits a continuum limit, 
the UV behavior of its Fourier transform Gyad'^) needs 
to be further analyzed. In particular, dyGyad'^)\T=o con¬ 
tains divergences. One is a trivial contact term which 
appears in the relation d'^Gyad'^)\T=o = —5{t = 0) -b 
dtdt'Fp°‘'^{t,t')\t=t'■ After this contact term has been 
subtracted, there remains a divergence in 



where D denotes the correlator of the field x and with ini¬ 
tial correlations = X],, Pl^Pi^P^^addxv^d^P-F -p). 


( 9 ) 
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where we used Eq. Q. This leading behavior at large A 
is not modified at non-zero T,p or m^. This correlator 
contributes to the energy density |15j . and in equilib¬ 
rium the divergence can be removed by a cosmological 
constant counterterm. 


RENORMALIZATION OUT OF EQUILIBRIUM 

To study the behaviour of Fp{t, t') for A —)• oo, we note 
that a formal analytical solution for Fp{t,t') is given by 
Fp{t,t') + with 


hom 


F, 


rriinh.G 


= Ppit)Ppit')dtdt'Fp{0, 0) -t ap{t)ap{t')Fp{0, 0) 
+ {(^p{t)pp(t') + ap{t')pp{t)) dtFp{0, 0), 

= I du I dv Pp{t,u)Ilp^p{u,v)pp{v,t'), 


/o 


/o 


irinh.nG ^ 

Pp = 4 


du pp{t,u)Ii^‘i{u)pp{Q,t') 


'0 




dupp(t,0)n;“(i;)pp(z;,t'). 


( 10 ) 


where we introduced Pp(t) = pp(t,Q) as well as <Jp{t) = 
—dt'Pp{t,0) + J* Pp{t,u)ll^p{u). For a Gaussian ini¬ 
tial state, vanishes identically and ap{t)^°‘'^^^ = 

—dt'Pp{t,0), in accordance with [15]. 

We first investigate the contribution in 


(101 which is independent of the initial conditions. Po¬ 


tential UV divergences can arise only from the vacuum 
part of Hi? (i.e. n,, Up-q —>■ 0) because the thermal con¬ 
tribution is exponentially suppressed for large loop mo¬ 
menta. Keeping only this part and using the Fourier 
representation Pp{t) = J one obtains 


T^inh,G p\ _ ^ Vo f ^ ^ 


(e*W* - e*‘^‘)(e 


dw 

^ J 

— iVad 


doj' 

2tt 


Pp{uj)pp{uj') 


') 


, (11) 


4cUqUJp—q(f}q CL?)(Qg UJ ) 

where flq = ujq + ojp-q- Note that the integrand has no 
poles because the numerator vanishes for —>■ uj,aj', 
respectively. The integration over q is superficially loga¬ 
rithmically divergent. To extract the most UV sensitive 
terms we use / duipp{uj) = 0 to rewrite the integral in an 
equivalent form, with the second line of 0 replaced by 




AlOqUJp—q^q 


1-h 


flq{uj + Ul') — UJU}' 

{flq - Lu){flq - UJ') 


( 12 ) 


Using Pp{u}) oc oj~^ for large w, one shows that the 
second term in the square bracket of (12) leads to ab¬ 
solutely convergent contributions to dtFp^^F 


and dtdt'Fp^^F^ Potential divergences therefore can 
only arise from the first term in this bracket. Using 
/ ~ we obtain 


F;^^^^{t,t') = f Pp{t)pp{t')Lp- pp{t)Sp{t') 


1 

TT L' 


~Pp{t )Sp{t) Cp{t — t) -l-... (13) 

where the ellipsis stands for absolutely convergent con¬ 
tributions, and we defined the integrals 


= 327r 


Sp{t) = 32 tt^ 


Cp{t-t') = 327r2 


d^^Q 


1 


^q^p—qf^'q 


{2tt) 

d'^q sm{flqt) 
^q^p—q^^ 


In — 
rriq, 


{277)° 

d^q cos{flq{t — t')) 
{2Tr)° 


^q^p—qf^q 


(14) 


For d = 4, Lp is logarithmically divergent for large A, 
while Sp and Cp are absolutely convergent for all t,t'. 
Nevertheless, 5p(0) = Lp and C'p(O) = —Lp are loga¬ 
rithmically divergent, which affects the correlators dtFp 
and dtdt'Fp (see below). The term oc Cp in dtdfFp'^^F 
matches the logarithmic divergence of the correspond¬ 
ing vacuum correlator for equal times ([^. In the fol¬ 
lowing, we discuss how these divergences affect the non¬ 
equilibrium correlators and demonstrate explicitly how 
they can be removed by the homogeneous and non- 


Gaussian contributions in (10) for a proper choice of ini¬ 


tial conditions. Before that, we briefly discuss the Gaus¬ 
sian case. 


Gaussian initial condition 

On general grounds, one expects that a physical ini¬ 
tial state should differ from the vacuum correlations by 
a finite, cutoff-independent amount. Implementing this 
idea rigorously would require to take initial n-point cor¬ 
relations into account for all n. In practice, one has to 
cut at some finite n. Let us first consider the Gaussian 
case 


(Gl) 


Fp(0,0) 

5tFp(0,0) 

dtduFp{0,0) 


Tpvac 




(0,0) + A(o) 

( 0 , 0 ) 


Ad) 

p 


= at9t-F;“(0,0) + 
= 0 for n > 3, 


Ad) 


(15) 


where only the connected two-point function is non-zero 
initially, and Ap) are cutoff-independent functions that 
parameterize the non-equilibrium initial state. The log¬ 
arithmic divergence contained in 9(5(/Fp“'=(0,0), cf. (|^, 
leads to a logarithmic divergence in the homogeneous so¬ 


lution (101, 


F!^<^'^{t,t') = ^-pp{t)pp{t')\n — 
TT m. 


- p 


■ finite. 


(16) 




















4 


This divergence has precisely th e sa me time-dependence 
as the one oc Lp in , cf. (13), but when summing 


both contributions there is in fact no cancellation. In¬ 
stead both divergences add up, and therefore the choice 
(Gl) does not admit a continuum limit for Fp{t^ t'). This 
can also be seen in the numerical solution, shown in Fig.[^ 
(dashed lines). 

Is it possible to remedy this shortcoming without going 
beyond the Gaussian initial state? To answer this ques¬ 
tion, we consider an alternative initial condition for the 
mixed derivative (and with the other derivatives initial¬ 
ized as in (Gl)) where we add ‘by hand’ a piece that re¬ 
moves the logarithmic divergence in Fp{t, t') at all times. 


(G2) dtdt,Fp{Q,Q) = 9tat,F™=(0,0)-2^Lp+A(2). ( 17 ) 


Indeed, Fp(t,t') possesses a continuum limit, as can also 
be observed in Fig.(upper graph, dotted lines). How¬ 
ever, closer inspection shows that this choice leads to the 
appearance of initial-time singularities in the two-point 
functions involving the canonical momentum, in partic¬ 
ular from ( |I^ and 0 it follows that 

dtF^°^\t,t')\t^o = --Pp{t')Lp + ... (18) 

TT 

is logarithmically divergent for A —>■ 00, as our numeri¬ 
cal simulation also confirms (not shown). Moreover, the 
momentum-momentum correlator 


dtdvFf^\t,t') 


.1 
TT L’ 


Pp{t')Sp{t) + pp{t)Sp(t') -\- Cp{t 


f —^Lp -|-... t,t' ^ 0 
(-l-^Tp-|-... t — y ^ If A ’ 



exhibits a cutoff-dependent ‘jump’ from the initial value 
imposed at t = t' = 0 to the value that matches the 
vacuum correlator (|^ at ‘late’ times (see Fig.[^ lower 
graph, dotted lines). Note that both (Gl) and (G2) lead 
to a cutoff-dependence which cannot be removed by a 
cosmological constant counterterm in the contribution of 
the field- and momentum correlator to the energy density, 
respectively. 


Non-Gaussian initial condition 




FIG. 2: Time-evolution of the statistical propagator Fp{t,t) 
and the momentum-momentum correlator dtdt'Fp{t,t') — 
dtdt'Fp°'‘^{t,t')\t^ti (for p = 0) at equal times for three differ¬ 
ent values of the cutoff Afm^ = 10, 25,100, and three differ¬ 
ent initial conditions (Gl) (dashed), (G2) (dotted) and (NG) 
(solid). We used A^q = = OiAp^g = n‘”m^ 

with n*" = 0.5, 'yfm^ = 0.28, = 0. 


avoids the pathologies in the Gaussian case and admits a 
well-behaved continuum limit. Using the matching pro¬ 
cedure developed in [14], 


In the following we demonstrate that the non-Gaussian 
initial condition 

(NG) Fp(0,0) = f;“^(0,0) + AW 

atFp(0,0) = 5tF;“=(0,0)-hAW, 
dtdt>Fp{0,0) = dtdt>F;-%o,o) + A^^\ 

03 = 03 “'”, On = 0 for n > 4, (20) 

characterized by two-point functions as for (Gl) and an 
initial three-point correlation equal to the one in vacuum 


la 


ijk^vac 


-2A 


0J„ + OJp-n + Vp 


( 21 ) 


for ijk = sss,aas, saa. All higher n-point functions are 
set to zero initially. 

The inhomogeneous part of the solution (10) now con¬ 
tains an additional piece involving a^. An analogous 
computation as above shows that 


= l\-2ppit)ppit')L 

TT L 


( 22 ) 
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~\~Pp{t)Sp{t') + pp{t')Sp(t) + ... 


Remarkably, the term cx Lp has the same structure as in 
t'), c/. (131, but with a relative factor —2. To¬ 


gether with the divergence in the inhomogeneous Gaus¬ 


sian part (13), this is precisely what is needed to cancel 


the logarithmic divergence of the homogeneous part ( 161 . 


In addition, it is important to note that the terms pro¬ 
portional to Sp cancel with those in 

This has several consequences which we want to stress: 
(i) Fp{t,t') and dtFp{t,t') converge to a finite contin¬ 
uum limit, (ii) dtdt'Fp{t^t') has a time-independent log¬ 
arithmic divergence for t = t' which matches precisely 
the one in vacuum, i.e. the difference dtdt'Fp{t,t') — 
dtdt'Fp°'‘^{t,f) also converges for all t,P > 0 . (hi) there 
are no initial-time singularities. These features can be 
observed also for the numerical solutions (see Fig.[^ solid 
lines), which are almost indistinguishable when varying 
the cutoff. Furthermore, (i) and (ii) imply that the en¬ 
ergy density is finite at all times and renormalized by 
the same counterterms as in equilibrium. We emphasize 
that the initial three-point correlation sizeably affects the 
solution Fp{t,t') not only at early times, but up to the 
thermalization time-scale t ~ I/ 7 . 


CONCLUSION 


We have shown for the first time how a proper account 
of initial, non-Gaussian vacuum correlations yields 
UV finite time evolution of the field- and momentum 
two-point correlator for all times, starting from an initial 
state that can be arbitrarily far from equilibrium, within 
a non-secular approximation scheme that captures 
thermalization at late times. The scheme is based on an 
expansion of initial correlations around the interacting 
vacuum state, complemented by usual vacuum coun¬ 
terterms. It is well-suited for analytical and numerical 


evaluation and allows a straightforward generalization 
to more complex theories, opening the way to the for¬ 
mulation of a renormalized initial value problem in QFT. 
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